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ABSTRACT 

Determining the preferred spatial location of the energy input to solar coronal loops 
would be an important step forward towa rds a more comp l ete u nderstanding of the 



coronal heating problem. Following on from lSarkar &: WalshI (|2008l ) this paper presents 



a short (10^cm= lOMm) "global loop" as 125 individual strands, where each strand is 
modelled independently by a one-dimensional hydrodynamic simulation. The strands 
undergo small-scale episodic heating and are coupled together through the frequency 
distribution of the total energy input to the loop which follows a power law distribution 
with index ~ 2.29. The spatial preference of the swarm of heating events from apex to 
footpoint is investigated. From a theoretical perspective, the resulting emission mea- 
sure weighted temperature profiles along these two extreme cases does demonstrate a 
possible observable difference. Subsequently, the simulated output is folded through the 
TRACE instrument response functions and a re-derivation of the temperature using dif- 
ferent filter-ratio techniques is performed. Given the multi-thermal scenario created by 
this many strand loop model, a broad differential emission measure results; the subse- 
quent double and triple filter ratios are very similar to those obtained from observations. 
However, any potential observational signature to differentiate between apex and foot- 
point dominant heating is possibly below instrumental thresholds. The consequences of 
using a broadband instrument like TRACE and Hinode-XRT in this way are discussed. 

Subject headings: Sun - activity; Sun - corona; Sun - UV Radiation; hydrodynamics; 
methods - data analysis 



1. Introduction 

Considerable effort has been devoted to date to determine the nature of solar magnetic plasma 
loops, the building blocks of the solar atmosphere. Particular interest has centered upon the 
examining of any temperature variation along an observed loop which could be related, through 
model comparison, to a possible prescribed heat input that has a well defined spatial localisation 
associated with it. Thus, there are two important avenues of investigation linked together here. 
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Firstly, one can take an observed t hermal profile say and en deavour to reproduce the temperature 
structure via a hydrostatic model. iPriest et al.l (|l998l . 12000 ) outlines evidence that the heat input 
in their large (> 700Mm) Yohkoh soft-Xray loop wa s distributed uniformly in the upper corona l 
portion (1.6 - 2.2 MK). This has been challenged bv Uschwanden et al.l (|200ll ) and|Reald (12002), 
who argued for footpoint and apex dominant heating respectively from the same dataset. This has 
been further investigated by lAdamakis et al.l ([2003) who, using Bayesian analysis methods, also 
argues fo r footpoint heating. Of parti cular relevance to the work outlined in this paper are the 
results of IWinebarger WarrenI (j200J), who studied the dynamic behaviour of the temperature 
and density of an impulsively heated ID hydrodynamic loop. From their analysis they conclude 
that TRACE observations of cooling loops are unable to provide ade quate informa t ion a bout the 
magnitude, duration and location of the earlier energy release. Also, iParenti et al.l (j2006l ) employ 
forward modeling techniques upon a nanoflare-heated multistranded loop to suggest that only high 
temperature coronal lines are capable of "remembering" the original heating signature. 

Secondly, one can take a theoretical model and fold the resulting plasma parameters through 
the instrument response functions to mimic what a given telescope would observe. With that 
in mind, it is very important to be aware of what current observational results are showing us. 
Consider concentrating upon the in t erpretation of EUV ima g ing data only . A nu m ber of authors 



( Aschwanden Nightingale! (|2005l ^. iNoglik fc WalshI ([20071), iNoglik et al.l ([2008|), ISchmelz et al 



20071 )) have used TRACE filter ratio data in an attempt to investigate t he physical pr o pertie s 



of coronal loops wi th mixed success . Whe t her one uses d ouble filter (e.g. ISchmelz et al.l ([2003)) 
or triple filter (e.g. iNoghk &: WalshI (j2007l ). INoglik et al.l (120081 )) ratios, it appears to be the case 
that a consistent, unique temperature determination is not possible. On top of that, another 
important factor is whether a loop is a single, monolithic structure or is the amalgamation of many 
sub-resolution evolving strands. It has been observed that the cooling time for an EUV loop is 
longer than the cooling time for a monolithic plasma structure (e.g. IWinebarger et al.l (|2003l ). 
Ugarte-Urra et al.l (|2009l )). hence leading to the suggestion that this extended cooling period is 
the result of a. bundle of successive, locally impulsively heated (and then cooling) plasma strands 
(|Warren et al.1 ^20od )). 



The concept of a l oop a s consisting of numerous su b-resolution elements was investigated ini- 
tially by Cargilj ( 1994 ) and Cargill &: Klimchuk ( 1997 ). where the authors assumed each strand 
could be represented by a single temperature and density value only. By heating each strand with 
several nanoflare-like energy events and allowing them to cool down by conduction and radiation, 
these au thors constructed and inyestiga. ted the behaviour of a multi-stranded global loop. Subse- 
quently, iPatsourakos &: KlimchukI (|2005l ) used their nanoflare heated multistranded model to show 
that the TRACE an d Yokoh S XT e mission was only weakly affected by the various dominant heat 
deposition locations. IWarrenI (|2006l ) also used the nanoflare heated multi-stranded model to show 
flare lightcurve can be better explained when the loop is divided into subresolved strands and 
heated impulsively. 



In this paper, we consider once again the multistrand loop model introduced in lSarkar &: Walsh 
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(|2008l ) (heareafter SW08) but here investigate the physical and observable consequences of spatially 
localising the nanoflare-heat input to either the loop structure apex or footpoint areas. Section 
2 outlines briefly the multistrand model while in Section 3, the results for apex and footpoint 
dominant heating are outlined. The observed consequences of such modelling are describ ed in 
Section 4, c o ncent rating upon recreating the filter ratio methods from lNoglik &: WalshI (j2007l ) and 



Reale et al.l (120071 ). while further work in this area is outlined in Section 5. 



2. Numerical simulation of multi-strand loop 

2.1. Model description 

The particular nanofiare multi-strand model employed in this paper was developed and outlined 
in SW08. Here we describe briefly the important aspects. The model simulates a short loop of length 
10 Mm. This loop consists of 125 strands where each individual element evolves hydrodynamically 
independently from each other; that is, there are 125 separate one dimensional hydrodynamic 
simulations. However the strands are coupled together through the frequency distribution of the 
total energy inp ut across th e loop s. Our ID hydrodynamical code is based on the Lagrange- 



Remap method (lArber et al.l (j200ll )). which is able to capture the shock fronts developed during 



the episodic heating events. 

At the beginning of every simulation all strands are at 10^ K and the plasma velocity is zero. 
The pressure and density maintains an exponential proflle to represent gravitationally stratifled 
plasma. There exists a highly dense plasma region (~ IMm long) at each end to mimic the 
chromosphere in the usual way. 



2.2. Heat input to the strands 



It is widely accepted that the often subtle interactions between the solar atmospheric magnetic 
fleld and the photospheric motions into which it is embedded, play the vital roles in determining 
a potential coronal heating mechanism. The question arises as to which of the compe ting physical 
timescales (driving versus alfven) is dominant and hence whether AC (wave heatin g: 



Muller et al 



(|1994| )) or, DC (magnetic reconnection of a tangled fleld or nanoflares: iParkeii (|l983l . 119881 )) heating 
is taking place. In this paper we concentrate upon the DC mechanism of small, localised energy 
bursts occurring along a loop structure. 

Hudson (1991) suggests that in order for nanoflares to provide the dominant heat input to the 
corona, the power law index associated with the frequency distribution of the energy event size 
should be more than 2. In SW08, we studied the effect of varying this index from to 3.29. In 
the following simulations we keep a = 2.29 (Case B in SW08) with a nanoflare energy range of 
10^^ to 5 X 10^^ erg. Figure H^a) displays a histogram of event size versus occurrence rate for the 
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simulations. In SW08, the event locations were spread randomly along the strands between —3.9 to 
3.9Mm- this produces an effective spatially uniform heat input to the loop. In these investigations, 
we now localise spatially where the events can occur. 

We consider two extreme cases as shown in Figure[T]^b) which is a "location" histogram. Firstly, 
there is the "footpoint dominant" heating (FDH) case where the vast majority of the nanoflares 
occur within a narrow region of ±3.5 to ib3.9Mm. Secondly, there is the case where the episodic 
energy events are clustered close to the loop apex at s = (termed "apex dominant" heating or 
ADH). Note that (i) in all cases, the overall total amount of heat provided to the strands is the 
same- the main investigation here is the impact on what we could observe when the heat input is 
distributed differently; (ii) the nanoflare lifetimes are chosen randomly over a range of 50 to 150 s 
and release their energy over a length scale of 0.02 Mm. The consequences of this spatial variation 
in the heat input are discussed in the following sections. 

3. Results of multi-strand model 
3.1. Effect on loop apex temperature 

As a reminder, SW08 defines the emission measure weighted temperature (Tem) as 

125 

^n^{s,t)dl{s)Ti{s,t) 



i=l 

where n indicates particle density, dl is the physical length of the computational grid, T is 
the temperature, s the position along the loop, t is the time and i is the number of the strands. 
Subsequently, Figure [2] plots Tem at the loop apex for three distinct cases. Note that Figure [2]^b) 
displays the same plot as Figure 6 (top left) in SW08 which is the case of spatially uniform heating 
(UNI). Of specific interest in SW08 the sudden dips in temperature, shown to be related directly to 
cool plasma blobs travelling along individual strands. It was demonstrated that these are generated 
when nanoflares release their energy further down the legs of the strands towards the chromospheric 
region. This in turn causes chromospherically evaporated blobs which travel along the strand and 
hence effect the calculation of Tem- 

Figure shows the ADH case. There are two aspects to note. Firstly given the spatial 
preference for the nanoflare bursts to occur far from the loop legs, no plasma blobs are created 
along any of the 125 strands. Secondly, there is an increase in the apex Tem up to ~ 2.4MK 
compared to ~ 2.2MK in the uniform heating case. In contrast. Figure [2jc) displays the FDH case. 
Compared to the UNI case, there are more plasma blob events influencing Tem^ should be the 
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case given the preferred energy deposition locations. Also, the average Tem is depressed further to 
~ 2.1MK. Examining the resulting Tem differences further, Figure [3] plots the time-averaged Tem 
profile along half of the loop length [-5, 0]Mm for ADH and FDH over a period of ~ 1500s from 
a start time of 7.3 x lO^s- this time window is chosen to avoid any major plasmoidal flows along 
the strands. We notice that there is a distinct diff erence in these thermal profiles. As mentioned 
by other authors (jPriest et al.l (|200d ) IWalshl (|l999l )). the preferred location of the dominant heat 
input to the strands (and subsequently to the amalgamated loop) does have a measurable effect on 
the averaged resulting temperature structure. 

For the ADH case, the nanoflare heating is taking place in an environment where thermal 
conduction dominates as the energy mechanism. Subsequently, the average temperature and the 
temperature gradient will increase to "remove" the deposited heat. In contrast, for the FDH case, 
the thermal profile is much fiatter in the "coronal part" of the loop. The nanoflare energy deposition 
is taking place close to the strand/ loop legs where radiation will play a more important role as 
well as a limited conductive loss. 



4. EUV imaging observations and the rederivation of the loop temperature from 

filter-ratio method 

In recent years, it has become possible to undertake forward modelling of simulation data 
by folding a calculated temperature/density through the response function for our instrument of 
choice. Thus, let us consider producing the synthetic emission for three passbands on the Transition 
Region and Coronal Explorer (TRACE) mission- 171, 195 and 284 A. Initially 

h,i = Gx{TWi{s,t)ds (2) 

where is the intensity of the strand i at passband A. Gx{T) is the temperature response function 
at passband A, rii is the density of the strand i, which is a function of space(s) and time(t) and ds is 
the computational grid length. While calculating the intensity, Feldman (1992) coronal abundances 
and Arnaud & Raymond (1992) ionization balance are considered. After calculating the emission 
from all the strands, one can derive the overall intensity of the global loop as, 

125 

h = Y.h, (3) 

i=l 

Once the loop intensity for each passband is derived, we can employ different ratio techniques 
to see the effect of a multistranded model. However it must be kept in mind that one of the 
main assumptions behind all of the following ratio methods is the notion that the volume of plasma 
under investigation in an instrumental pixel is isothermal. This is not the case for our multi-strand, 
multi-thermal model. However, it is instructive to examine what the subsequent consequences of 
this in relation to how current observational datasets are being analysed. 
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4.1. Single filter ratio 



The single filter ratio method has been used by several authors (e.g. lAschwanden et al.l ([20001)) 
to estimate the temperature of a target loop from EUV imaging data. Using the synthesised EUV 
emission produced from our model as outlined above, Figure S] displays the time evolution of three 
filter ratios for the (a) ADH and (b) FDH cases for a TRACE-sized pixel (1 arcsec) placed at the 
loop apex. A time integration (exposure time) of 60s has been employed. 

There are several aspects to note. Firstly, compared to the Tem, none of the three filter 
ratios provided a "correct" estimate of the loop temperature. Whereas the average Tem for the 
ADH (FDH) case was ~ 2.4MK (2.0MK), the three filter ratios vary around ~ 1.8MK (284/195), 
~ 1.6 MK (284/171) and ~L25MK (195/171) respectively. Observationally it is always seen 



(e.g. (jSchmelz et al.ll2003l . 120071 )) from TRACE and FIT active region data that the temperature 
obtained from the 195 to 171 ratio is always statistically indistinguishable from 1.2MK and the 195 
to 284 ratio always biased towards 1.8MK. There appears to be no distinct observational difference 
in the filter ratio temperature between ADH and FDH, apart from the fact that the ADH variation 
appears to be larger, a natural consequence as the dynamic heating is taking place directly at the 
apex. 

Now consider deriving the thermal profile along the loop from these synthesised TRACE emis- 
sions. Employing the 284/195 A ratio with a 60s integration time. Figure [5] displays the resulting 
temperature variation along half the loop and slightly above the chromospheric boundary. Com- 
paring this to Figure EJ the resulting thermal profiles for ADH and FDH are different. However, 
the temperature discrimination of TRACE would not be adequate to detect this slight difference 
of ~ 0.1 MK, which would be of the order of the observed error bar. 



4.2. Color-Color method 



Given the apparent problems with t he single filter ratio , a double filter ratio me thod has been 
employed that combines all three ratios (jChae et al.l (|2002l ). iNoglik &: WalshI (|2007l )). Here single 
filter ratios are plotted against one another to create a unique temperature-dependent color-color 
(C-C) curve. In this case, we define 
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where Ci and C2 are now functions of temperature. Figure [6] plots the theoretical C1/C2 curve with 
the corresponding values derived from (a) ADH and (b) FDH scenarios, where the same TRACE 
pixel-sized portion of the loop apex is used over an integration time of ~ 1500s. In this case, a 
much longer integration time of ~ 1500s is used to include any longer term variability. 
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In both cases, the synthesised data points cluster in a region far from the C-C curve. As the 
bounding box shows, there is a sUght difference between ADH and FDH, corresponding to a smah 
decrease in both ratios. This is due to the perceived small decrease in "temperature" indicated in 
Figure O 

It should be noted that this type of cl ustering behaviour has been observed in TRACE obser- 
vational data. For example, in Figure 5 of lNoglik &: WalshI (j2007l ). there is a clear demonstration 
of clustering of datapoints along loop structures observed at the solar limb. 



4.3. Combined filter ratio method 

Given the above unsuccessful attempts to predict a temperature from single and double fil- 
ter ratio ni e thods , it is instructive to consider a combined filter ratio (CFR) as introduced by 



Reale et al.1 ((20071). 



It must be noted that those authors employ this method upon XRT observations where they 
have many more filters available compared to the three TRACE passbands considered here. Briefly, 
CFR is the ratio between the geometric mean of the emission detected in all available filters and 
the emission Ix in a single (A-th) filter given by 



(11 ^a)^/" (n^Am)^ 



/n 



CFR,(T) = = ^--g^,^) (6) 

Thus in this case, n = 3. A further improvement is obtained by taking the product of the 
ratios computed for the two filters. The resulting combined improved filter ratio is defined as 

CIFR{T) = CFRi{T) x CFR2{T) (7) 

Figure [T| plots the CIFR for the FDH only- there is no significant difference between these curves 
and the ADH case. Note that while plotting we again consider the average value over a TRACE- 
sized pixel placed at the loop apex (s = 0) location. The fluctuations observed during the first 
1000s are due to the high fiuctuations in density at the beginning of the simulation. CIFR predicts 
values for the temperature that are quite different than the values predicted by the corresponding 
filter-ratio temperature (Figure [Ub)). In this case, the calculated temperatures are about 1.65MK, 
1.5MK and 2.1MK for CIFR(171 x 195), CIFR(195 x 284) and CIFR(284 x 171) respectively. The 
fluctuations in the CIFR temperature evolution are less than the single filter-ratio as well. 
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5. Summary and discussion 

In this paper we have examined a multi-stranded loop model where the episodic heat input to 
the structure is deposited at two extreme locations; namely, at the strand/loop apex (ADH) and 
footpoint (FDH). The time averaged EM weighted temperature along the loop that results from 
the overall effect of all 125 strands was calculated. As shown in Figure [3l the ADH thermal profile 
has a steeper temperature gradient in the coronal (> 10^ K) part of the loop than that resulting 
from FD H. This mirro r s similar results fro r n hyd rostatic modelling of spatially dependent heating 
in loops dWalshl (|l99d ). kschwanden et al.l wi))- 



Subsequently, the multi-strand simulation values were folded through the instrument response 
function for three EUV TRACE passbands (171, 195, 284 A) and the cumulative intensity across 
all strands calculated. The resulting emission was then used to calculate a loop plasma temperature 
by three filter ratio techniques- single, double and combined improved. The resulting values and 
evolution of a TRACE-sized pixel at the loop apex is very representative of what we evaluate from 
TRACE loop observations. 

Also it is shown that, using the single filter ratio, the derived temperature along the length 
of the loo p for FDH and ADH will not be differentiable observationally. This conclusion is similar 
to that in iPatsourakos &: KlimchukI (j2005l ) who, using their loop model showed that TRACE &: 
Yohkoh-SXT observation of loops de p end o nly weakly upon the spa tial distribution of nanoflare 
heating. Also IWinebarger &: WarrenI (120041 ) and IParenti et al.l ([200a) point out the TRACE and 
low temperature spectroscopic observations (where radiation is the dominant cooling phenomenon) 
are unable to retain information about the nature and location of the original heat input. 

Considering each in turn, we see firstly that the single filter ratio technique gives three quite 
different tempe rature values from the three ratios (Figure H]). It must be noted that this has been 
investigated by lWeber et al.l (120051 ) who employed a flat differential emission measure distribution 
to mimic an "extreme" multithermal plasma. The paper found clearly that the single filter ratio 
method is very biased towards the ratio of the integrals of the temperature response function 
under these conditions. Essentially the same scenario is happening in the current simulations. 
Multiple strands are heating up and cooling down simultaneously such that the spread of the true 
temperature and density values generated from all the strands will be quite broad. Although it 
will not be completely flat, the resulting DEM stretches over a several decades of temperature (i.e. 
~ 10^ — lO^K); thus the single filter ratio method will tend towards isothermal values for each given 
line pair. 

Similarly, a complementary argument can be made concerning the double filter ratio method. 
If there is a tendency towards a certain value of the ratio from a single filter ratio, this will be 
reflected in a clustering of the points in the C-C curve, as demonstrated in Figure El 



Finally consider the CIFR approach. Employing a flat DEM model as a proxy for a multither- 
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mal plasma (as in lWeber et alj (|2005l )) we can write from 



^171 

CFR,,, (T) = (-^in X /195 X /284)^/^ ^g^^ 



Since, 

/A~y" Gx{T)DEM{T)dT = DEM flat J Gx{T)dT 
then following equation ([7]) we will have the CIFR for 171 and 195 as, 



CIFR(T) = (^m X X /284)V-^ 

A?! X /i95 

rT^b /-Tfc 



( / ' Gni{T)dT X / ' Gi95(r)dr x / ' G284(T)dr)2/3 (9) 



I ' Gni{T)dTx I ' GiQ^{T)dT 



Here, Ta and T}y are the starting temperature and end temperature for the integration. If Ta and 
Tft are fixed, CIFR(171 x 195) always maintains a fixed ratio i.e., 0.6, which corresponds to the 
temperature ~ 1.7 MK the value we obtain from the simulations. The same is true for all other 
combination of filters. Thus like the single filter ratio, CIFR also gives the ratio of areas under 
the temperature response function curves. Hence, one can conclude that for a broad multithermal 
plasma, CIFR also will provide a flat temperature profile, that does not reflect the real thermal 
structure. 

In this paper we have produced the CIFR using only three passbands of TRACE. Thus even 
though our present result is restricted to within the TRACE passbands, the theoretical calculation 
will be true for Hinode-XRT passbands as well. A detail calculation considering the XRT passbands 
is in progress for our future work. 

Overall, this analysis demonstrates that the emission and filter ratio results we are getting 
from imagers can be reconciled with a multi-strand, multi-thermal loop approach. However the 
temperature values they produce do not correspond at all to the EM weighted temperatures, but 
are purely an instrumental effect. 
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Fig. 1. — (a)Energy histogram for the size of a heating event versus the number of times it occurs 
during the simulation across all 125 strands. This has been fitted with a straight line, to show 
the power law slope of a = 2.29; (b)Location histogram displaying the preferential location of the 

nanoflare events within the loop. Two extremes are considered -apex ( ) and footpoint ( ) 

dominant heating. 
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Fig. 2. — Time evolution of Tem at the loop apex for three cases (a) apex dominant heating (ADH) 
(b) spatially uniform heating (UNI) and a (c) footpoint dominant heating (FDH). 
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Fig. 3. — Time averaged emission measure weighted temperature profile along the half loop length 
ADH ( ) and FDH ( ) over a time period of ~ 1500s. 
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Fig. 4. — Single filter ratio temperature evolution at the loop apex for (a) ADH and (b) FDH loop. 
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Fig. 5. — Temperature profile along the length of the loop derived using 284/195 filter ratio. 




Fig. 6. — Simulated data points in green on the C-C plot for (a) ADH and (b) FDH loop. Note 
that the solid line is the C-C temperature curve. 
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Fig. 7. — Combined improved filter ratio temperature evolution using three TRACE passbands for 
the FDH case. 
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careful reading of the manuscript. 
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